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Abstract. We present four large sets of evolutionary 
tracks for stars with initial chemical compositions [Y = 
0.250, Z = 0.008], [Y = 0.273, Z = 0.019], [Y = 
0.320, Z = 0.040] and [Y = 0.390, Z = 0.070] and en- 
hancement of a elements with respect to the solar pat- 
tern. The major improvement with respect to previous 
similar calculations is that we use consistent opacities 
- i.e. computed with the same chemical composition as 
adopted in the stellar models - over the whole relevant 
range of temperatures. For the same initial chemical com- 
positions [Y, Z] and otherwise identical input physics we 
present also new evolutionary sequences with solar-scaled 
mixtures of abundances. Based on these stellar models we 
calculate the corresponding sets of isochrones both in the 
Johnson-Cousins UBVRIJHK and HST/WFPC2 photo- 
metric systems. Furthermore, we derive integrated mag- 
nitudes, colours and mass-to- light ratios for ideal single 
stellar populations with total mass equal to 1 Mq . Fi- 
nally, the major changes in the tracks, isochrones, and 
integrated magnitudes and colours passing from solar- 
scaled to a-enhanced mixtures are briefly outlined. Re- 
trieval of the complete data set is possible via the www 
page http : //pleiadi .pd. astro . it. 
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1. Introduction 

There are indications both from observations and theoret- 
ical studies that stars in many astrophysical environments 
may have the relative abundances of elements synthesized 
by nuclear a-capture reactions (e.g. C, O, Ne, Mg, Si, 
Ti, etc) that are enhanced with respect to the solar ones. 
In brief, King (1994) found the mean [O/Fe] ratio for a 
handful of halo stars to be +0.52 dex. Similar results are 
found for globular cluster stars (e.g. [a/Fe] = 0.18 dex 
for 47 Tuc, Carney 1996) and galactic bulge K giants 
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(Mc William & Rich 1994, Rich & McWilliam 2000). Ac- 
cording to models of chemical evolution of the Galaxy, 
bulge stars should be a-enhanced to reproduce the metal- 
licity distribution (Matteucci & Brocato 1990, Matteucci 
et al. 1999). Indications that the abundance ratios of ele- 
ments in elliptical galaxies are not solar (e.g. [Mg/Fe] = 
0.3 — 0.7 dex) come from Worthey et al. (1992). Also syn- 
thetic metal line strengths (Weiss, Peletier & Matteucci 
1995) from high-metallicity models show that solar-scaled 
models do not fit the observed Mg2 (5180 A) and Fe (5270 
and 5335 A) line indices. 

Given these premises, it is of paramount importance 
to include the enhancement of a-elements in the popu- 
lation synthesis models designed to describe those galax- 
ies. In the past this has been done but using models that 
suffered from many limitations. Firstly, till recently there 
were no extended sets of opacities for a-enhanced mixtures 
covering the complete relevant range of stellar tempera- 
tures. Such opacity calculations are nowadays available. 
Secondly, it has often been assumed that the behaviour 
of a-enhanced isochrones can be easily mimicked by suit- 
ably scaling the metallicity [Fe/H] (Salaris et al. 1993). 
However, the validity of such a re-scaling has only been 
confirmed for low metallicities, and it holds only for some 
particular cases like all a-elements being enhanced by the 
same factor. Weiss et al. (1995) already pointed out the 
differences in the evolution due to a-enhancement at so- 
lar and super-solar metallicities. Recent work by Salaris & 
Weiss (1998) and VandenBerg et al. (2000) also indicates 
that the simple re-scaling of the [Fe/H] content becomes 
risky at high metallicities. Moreover, we should keep in 
mind that there is no a priori reason, from the point of 
view of both nucleosynthetic and chemical evolution the- 
ories, to expect that all a-elements are enhanced by the 
same factor. From observations of halo and bulge stars, the 
degree of enhancement of different a-elements is, in gen- 
eral, not constant with differences amounting to ±0.2 dex. 

Therefore in this paper we intend to re-analyze the im- 
pact of the enhancement of a-elements on the properties 
of metal-rich stellar populations (whose metallicity ranges 
from half-solar to three times-solar). The analysis is made 
for a particular choice of a-enhanced abundance ratios, 
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for which self-consistent opacities are available. We start 
calculating extended sets of stellar evolutionary tracks 
(Sect, g), with total metal contents Z = 0.008, 0.019, 
0.040 and 0.070 both for solar-scaled (same abundance 
ratios as in the solar mix) and a-enhanced patterns of 
abundances. Particular attention is paid to the descrip- 
tion of the adopted initial element abundances and sources 
of opacity. Based on those stellar models we derive large 
grids of isochrones and single stellar populations for which 
we present magnitudes, colours and mass-to- light ratios in 
many pass-bands of the Johnson and WFPC2 photomet- 
ric systems (Sect. g). Extensive tabulations of stellar mod- 
els, isochrones and photometric properties are provided in 
electronic format for the purpose of general use. Finally, 
a short summary is given in Sect. |J 

VandenBerg et al. (2000) recently presented extensive 
sets of evolutionary tracks for a-enhanced mixtures, lim- 
ited to the range of low-mass stars (0.5 < M/M© < 1.0). 
A constant enhancement of all a-elements has been as- 
sumed (also in their opacity tables). In the present work, 
instead, the choice of initial abundances has been guided 
by observational results. It is worth remarking that, within 
the present errors, none of both choices seems to be 
strongly preferred. 

2. Stellar tracks 

We calculate four grids of stellar models with initial 
masses from 0.15 to 20 M Q and initial chemical compo- 
sitions [Y = 0.250, Z = 0.008], [Y = 0.273, Z = 0.019], 
[Y = 0.320, Z = 0.040], [Y = 0.390, Z = 0.070]. The 
metallicities and the helium-to-metal enrichment ratio 
AY/AZ are chosen in such a way to secure consistency 
with the Girardi et al. (2000) models. The helium-to-metal 
enrichment law is Y = 0.23 + 2.25 Z. 

For the above values of [Y, Z] , we compute both solar- 
scaled and a-enhanced sets of stellar tracks. The evo- 
lutionary phases covered by the grids extend from the 
zero age main sequence (ZAMS) up to either the start of 
the thermally pulsing asymptotic giant branch (TP-AGB) 
phase or carbon ignition. 

2.1. Input physics 

The input physics of the present stellar models is the same 
as in Girardi et al. (2000), apart from differences in (i) 
the opacities, and (ii) the rates of energy loss by plasma 
neutrinos. A short description of the updated ingredients 
and the way we deal with the enhancement of a-elements 
is presented in the sections below. 

2.1.1. Initial chemical composition 

For a fixed total metallicity Z, the solar-scaled abundance 
ratios of metals are taken from Grevesse & Noels (1993). 
In the a-enhanced case, we adopt a mixture in which only 




l°g T eff 

Fig. 1. Evolution of a 1 M© star with the indicated chem- 
ical composition up to the age of the Sun according to dif- 
ferent opacity sources and enhancements of a-elements. 
See text for the meaning of the various labels and more 
details. 

the elements resulting from nuclear a-capture reactions 
(O, Ne, Mg, Si, S, Ca, Ti, Ni) are enhanced with respect 
to the solar abundance ratios. This is the same mixture as 
used by Salaris et al. (1997) and Salaris & Weiss (1998) to 
study the age of the oldest metal-poor and disk metal-rich 
globular clusters. For the sake of clarity, Table [i] lists the 
abundances for the solar-scaled and a-enhanced mixtures. 
Columns (2) and (4) show the abundance A c \ of the ele- 
ments in logarithmic scale, A c \ = log N c \/Nn + 12, where 
iV c i is the abundance by number. Columns (3) and (5) 
display the abundance by mass fraction X e \, relative to 
the metallicity Z. Columns (6) and (7) list the enhance- 
ment ratio (relative either to Fe or the total metal frac- 
tion) in spectroscopic notation. The enhancement factors 
come from the determinations of chemical abundances in 
metal-poor field stars by Ryan et al. (1991, cf. also Salaris 
& Weiss 1998). To derive the initial abundances of differ- 
ent isotopes we use the isotopic ratios (by mass) of Anders 
& Grevesse (1989). 

2.1.2. Opacities 

The opacities we have adopted are the same as in Salaris, 
Degl'Innocenti & Weiss (1997). They are available at: 
http : //www-phys . llnl . gov/V_Div/0PAL/exist ing . html 
Details are given below. 

For high temperatures (log T > 4) we use the radiative 
OPAL opacities by Rogers & Iglesias (1995) in the case of 
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O 


8.870 


0.482273 


9.370 


0.672836 


0.50 


+0.1442 


Ne 


8.080 


0.098668 


8.370 


0.084869 


0.29 


-0.0658 


Mg 


7.580 


0.037573 


7.980 


0.041639 


0.40 


+0.0441 


Si 


7.550 


0.040520 


7.850 


0.035669 


0.30 


-0.0558 


S 


7.210 


0.021142 


7.540 


0.019942 


0.33 


-0.0258 


Ca 


6.360 


0.003734 


6.860 


0.005209 


0.50 


+0.1441 


Ti 


5.020 


0.000211 


5.650 


0.000387 


0.63 


-0.1875 


Ni 


6.250 


0.004459 


6.270 


0.002056 


0.02 


-0.3371 


C 


8.550 


0.173285 


8.550 


0.076451 


0.00 


-0.3553 


N 


7.970 


0.053152 


7.970 


0.023450 


0.00 


-0.3553 


Na 


6.330 


0.001999 


6.330 


0.000882 


0.00 


-0.3553 


Cr 


5.670 


0.001005 


5.670 


0.000443 


0.00 


-0.3557 


Fe 


7.500 


0.071794 


7.500 


0.031675 


0.00 


-0.3557 



solar-scaled mixtures, and by Iglesias & Rogers (1996) for 
the a-enhanced ones. At low temperatures these are com- 
bined with the molecular opacities provided by Alexander 
& Ferguson (1994) for the solar mixtures, and by Alexan- 
der (priv. communication) for the a-enhanced ones. In 
both cases, low- and high-temperature opacity grids are 
merged by interpolation in the temperature range from 
logT = 3.8 to logT = 4. For very high temperatures, 
logT > 8.7, the Weiss et al. (1990) opacities are used. 
These rely on the Los Alamos opacity database. The ra- 
diative opacities for C, O mixtures, needed in the post- 
main sequence phases, are from Iglesias & Rogers (1993). 
Finally, conductive opacities of electron-degenerate mat- 
ter are from Itoh et al. (1983). 

It is very important to emphasize that the present 
models are based on complete and consistent opacities, i.e. 
the opacities (both for solar and a-enhanced mixtures) are 
generated for the same pattern of abundances adopted in 
the stellar models, both at low and high temperatures. 

The importance of this issue is illustrated by Fig. pi 
which shows the effect in the Hertzsprung- Russell (HR) di- 
agram of adopting different (inconsistent) opacities and/or 
mixtures and/or enhancement factors as input to a 1 Mq 
main sequence model with [Z = 0.004, Y = 0.254]. Some 
of the illustrated cases are often found in literature. The 
tracks are evolved up to the solar age of 4.6 Gyr. The 
meaning of the labels is the following: the first label refers 
to high temperature opacity (logT > 4) and can vary 
from OPAL-S, i.e. solar-scaled mixture, to OPAL-a, i.e. 
enhanced according to Salaris & Weiss (1998; W95) or 
Rogers (priv. communication; R95). This latter case is 
a special a-enhanced mixture characterized, if compared 
with Salaris & Weiss (1998), by more C, less O, more Ne, 
more Na (which is not an alpha-element), less Al, more 
S, and more Fe. The second label refers to the molecu- 
lar opacity, with the same notation as for solar and a- 



enhanced mixtures. AF stands for Alexander & Ferguson 
(1994). The prescription enclosed in a box represents the 
one adopted by us. 

Another point to emphasize is the difference between 
the present opacities and those used by Girardi et al. 
(2000). Two are the main sources of differences. First the 
adoption of the Itoh et al. (1983) conductive opacities in- 
stead of the older ones by Hubbard & Lampe (1969). Inter- 
estingly, Catelan et al. (1996) show that Itoh et al. (1983) 
conductive opacities should not be used in RGB models, 
because they are valid only for a liquid phase, and not for 
the physical conditions existing in the interiors of RGB 
stars. However, we verified that the effect of using the Itoh 
et al. (1983) opacities instead of the Hubbard & Lampe 
(1969) is to increase the core mass at the helium flash by 
0.006 M Q (see also Castellani et al. 2000 and Cassisi et al. 
1998), which is a negligible effect. 

The second difference lies in the numerical technique 
used to interpolate within the grids of the opacity tables. 
In this paper we use the two-dimensional bi-rational cubic 
dampcd-splinc algorithm (see Schlattl & Weiss 1998; and 
Spath 1973), whereas Girardi et al. (2000) adopted the 
smooth bi-parabolic spline interpolation (SSI) already in- 
troduced in our code by Bressan et al. (1993). We verified 
that the differences amount to about a few percent. 

In general, evolutionary models of low mass differ very 
little by changing the opacity interpolation scheme. How- 
ever, massive stars (say M > 5Mq) are much more sen- 
sitive to the interpolation algorithm. This is shown in 
Fig. ||, where we plot two 4 M Q and 6 M© models with 
[Z = 0.008, Y = 0.250] calculated both with the SSI 
(continuous lines) and the damped-spline (dashed lines) 
scheme. 

Figure compares the effects of the different interpo- 
lations on the opacity at the physical conditions holding 
in the inner regions of a 6 Mq star when the central hy- 
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Fig. 2. Evolution of the 4 and 6 Mq models according to 
different opacity interpolation schemes. The dashed lines 
are tracks calculated with the damped spline, whereas the 
continuous lines are tracks calculated with SSI spline. The 
dots indicate the points when the central hydrogen mass 
fraction is X c = 0.2. See also the companion Fig. || for 
more details on the opacity. 

drogen mass fraction is X c = 0.2, i.e. at the stage marked 
with a dot in Fig. |^. At this stage the temperature at the 
Schwarzschild border of the convective core is logT = 7.38 
and the density is logp = 1.05 (corresponding to the ar- 
row in the Fig. ^). For the temperatures indicated along 
the curves, these represent the opacity as a function of the 
quantity logi? = logp — 3 logT + 18 (see Rogers & Igle- 
sias 1992 for more details). The continuous curve is based 
on the SSI interpolation, whereas the dashed curve refers 
to the method adopted in the present study. The small 
difference in opacity at the core border is typical of the 
values found throughout the evolutionary sequences, and 
is responsible for the different HR patterns. 

We notice that both interpolations are acceptable, be- 
cause the percentage error associated to both of them is 
roughly 3 %, which is lower than the opacity uncertainty 
in the original OPAL tables (the latter can be as high as 
10 % at individual T/p points). This shows that even the 
fine OPAL grid for some applications is not sufficiently 
dense and more grid points might be needed. This is also 
illustrated by the Solar Model by Schlattl & Weiss (1998), 
where the inverse situation, i.e. spline producing a nicer 
opacity run, is found. 

The reason for the different luminosities in the high- 
mass models, is that in these stars the radiative gradient 
across the border of the convective core is shallow. There- 



Fig. 3. The opacity k as a function of the density and tem- 
perature as indicated. The dashed lines are for damped- 
splined opacity interpolations, the continuous lines are 
for the bi-parabolic splined ones. The arrow shows the 
location of the Schwarzschild border of the convective 
core in a 6 Mq model with initial chemical composition 
[Z = 0.008, Y = 0.250] during the main sequence phase, 
when the central hydrogen mass fraction is X c — 0.2. 
For this model the temperature at the core border is 
logT = 7.38. 



fore, even slightly larger opacities at the Schwarzschild 
border may cause an increase of the convective core and 
lead to higher luminosities. This trend already begins dur- 
ing the core H-burning phase and remains during the sub- 
sequent evolution. The effects are detectable starting from 
stars with mass larger than 6 — 10 Mq (see Fig. |^) and 
becomes dramatic over 60 Mq. It is worth recalling that 
the maximum mass of our stellar models is 20 Mq . 

2.1.3. Neutrino losses 

Energy losses by neutrinos are from Haft et al. (1994). 
Compared with the previous ones of Munakata et al. 
(1985) used by Girardi et al. (2000), neutrino cooling dur- 
ing the RGB is more efficient. This causes an increase of 
0.005 Mq in the core mass at helium ignition. We have 
checked this on a 1 Mq model with Z = 0.019. Differ- 
ences of the same order have been found by other authors 
(e.g. Haft et al. 1994; Cassisi et al. 1998). 
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Fig. 4. Evolutionary tracks for the composition [Y = 
0.250, Z = 0.008] and the a-enhanced mixture of abun- 
dances. The initial mass in solar units is indicated along 
each curve. 



2.1.4. Remaining physical input 

The remaining physical input is the same as in Girardi et 
al. (2000), to whom we refer for further details. In sum- 
mary, we follow the evolution of H, 3 He, 4 He, 12 C, 13 C, 
14 N, 15 N, 16 0, 17 0, 18 0, 20 Ne, 22 Ne, Mg, according to the 
nuclear reaction rates of Caughlan & Fowler (1988) and 
Landre et al. (1990). For Mg we refer to the total content 
of 24 Mg, 25 Mg and 26 Mg. 

The equation of state (EOS) is that of a fully-ionized 
gas for temperatures higher that 10 7 K. For lower temper- 
atures we adopt the Mihalas et al. (1990 and references 
therein) EOS. 

Calibration of the solar model fixes the value of the 
mixing length parameter a of the Bohm-Vitense (1958) 
theory at 1.68. Overshooting both from the convective 
core and the convective envelope are taken into account ac- 
cording to the ballistic algorithm by Bressan et al. (1981). 
Table summarizes the values of the overshooting param- 
eters for the core, A c , and envelope, A we have adopted 
as a function of the initial mass. 

As far as the mass-loss is concerned, this is taken into 
account and its effect on the stellar structure is calculated 
for masses higher than 6 M Q . We used the same prescrip- 
tion adopted by Fagotto et al. (1994), where more details 
can be found. 



Fig. 5. Evolutionary tracks (same composition as Fig. [| 
for low-mass models up to the RGB-tip. 
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Fig. 6. Evolutionary tracks (same composition as Fig. |4|) 
for low-mass models from the ZAHB to the start of the 
TP-AGB phase. 



2.2. HR diagram 

For the sake of illustration, in Figs. |], || and ^| we 
present the complete set of evolutionary tracks for [Y = 
0.250, Z = 0.008] and a-enhanced mixtures. 
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Table 2. Adopted parameters in the overshooting model 
by Bressan et al. (1981) and Alongi et al. (1991), as a 
function of the initial mass M (in M & ). 



M/Mf 



© 



A e 



0.5 - 
1.1 - 
1.5 - 
2.0 - 
2.5 



1.0 
1.4 
2.0 
2.5 
20 




0.5(2.5 



0.25 
0.25 
0.25 
M) + 1.4(M 
0.7 



2) 



Fig. |i| shows the models of intermediate-mass stars 
from the zero age main sequence (ZAMS) up to the TP- 
AGB phase and those of massive stars from the ZAMS 
up to carbon ignition. Fig. || presents the low-mass tracks 
from the ZAMS up to the RGB-tip. Finally, Fig. | shows 
the corresponding He-burning phase from the zero age 
horizontal branch (ZAHB) up to the TP-AGB phase. 

Some relevant physical quantities are given in Table |, 
which lists, as a function of the initial mass, the lifetimes in 
years of the central H-burning and He-burning phases, in 
and £hc respectively, the core mass at the helium ignition, 
M c (HeF), and the core mass at the first thermal pulse 
on the AGB phase, M c (l st — TP). For intermediate-mass 
and massive stars, M c (HeF) refers to the core mass at the 
Hc-ignition. For the same chemical composition Table || 
reports the changes in the surface chemical abundances of 
those elements that suffered from nuclear processing via 
the pp-chain and CNO-cycle in the deep interiors, after 
the completion of the first and the second dredge-up by 
envelope mixing. The surface abundances of heavier ele- 
ments are not altered by nuclear processing. The last two 
columns give the ratios of 14 N and 16 with respect to 
their initial values. These ratios are indicative of the effi- 
ciency of the dredge-up episodes. 

Table [| displays the transition masses M conv , MhoF 
and M up for models calculated with solar and a- 
enhanced initial composition. In brief, M conv is the maxi- 
mum mass below which stars become fully convective dur- 
ing the central H-burning phase, A/hcF is the maximum 
mass for the core He-flash to occur, and finally M up is 
the maximum mass for the central C-burning to occur in 
highly electron degenerate gas (deflagration or detonation 
followed by disruption of the star). 

-^HcF is here defined as the initial mass for which the 
core mass at He ignition has its minimum value (see Ta- 
ble || and Fig. ^). Fig. ^ shows that similar values for MjjeF 
are found both in the a-enhanced and in the solar-scaled 
case. Compared with Girardi et al. (2000), we find higher 
Af up values. This is once more caused by the different in- 
terpolation schemes of the opacity tables. The one cur- 
rently in use yields smaller convective cores for the same 
initial mass. It follows that with the present overshooting 
prescription the maximum mass for a star to develop an 
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Fig. 7. Core mass at He ignition as a function of the 
initial mass, for both solar-scaled (upper panel) and a- 
enhanced (bottom panel) models. 

Table 5. The transition masses M conv , A/hcF and M up . 



Y 



mixture 



A I] 



HeF 



0.008 
0.019 
0.040 
0.070 



0.008 
0.019 
0.040 
0.070 



0.250 
0.273 
0.320 
0.390 



0.250 
0.273 
0.320 
0.390 



solar-scaled 
solar-scaled 
solar-scaled 
solar-scaled 



a-enhanced 
a-enhanced 
a-enhanced 
a-enhanced 



0.36 
0.37 
0.36 
0.33 



0.37 
0.37 
0.36 
0.33 



1.90 
2.10 
2.20 
1.90 



1.90 
2.10 
2.10 
2.00 



5- 6 

6- 7 

7- 8 
7-8 



electron-degenerate core after the core He-burning phase 
is slightly higher than previously estimated. 



3. Isochrones 

The four sets of stellar models are used to calculate 
isochrones, integrated magnitudes and colours of single 
stellar populations with ages from 10 7 to 2 x 10 10 yr both 
in the Johnson-Cousins UBVRIJHK and HST-WFPC2 
photometric systems. The age range is large enough to 
describe young clusters and associations as well as old 
globular clusters. The dense grid of stellar tracks allows 
us to construct detailed isochrones at small age steps 
(Alogt = 0.05). This also makes the data-base suitable 
for the simulation of synthetic colour-magnitude diagrams 
(e.g. Girardi 1999). 
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Table 3. Lifetimes of the central H-burning and He-burning ijje phases, the core mass at the helium ignition 
A/ c (HeF) and the core mass M C (1-TP) at first thermal pulse on AGB phase, for the models with Z — 0.008 and 
solar-scaled mixture (columns 2 to 5) and for the models with Z = 0.07 and a-enhanced mixture (columns 6 to 9). 
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= 0.008 solar-scaled 
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= 0.07 


a-enhanced 




M/Mq 


Wyr 
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2.45 10 10 
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l.UO 
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0.4434 
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0.4781 
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Before theoretical isochrones are constructed, the TP- 
AGB phase is included in all tracks of M < M up . There- 
fore, this phase is present in the isochrones older than 
about 10 8 yr. The reader is referred to Girardi & Bertclli 
(1998) and Girardi et al. (2000) for all details about the 
synthetic algorithm used to follow the TP- AGB evolution. 



3.1. Isochrones in the HR-diagram 

The HR diagram of Fig. || shows a comparison between 
isochrones with chemical composition [Z = 0.019, Y = 
0.273] both for solar-scaled (continuous lines) and a- 
enhanced (dotted lines) stellar models. For the sake of 
clarity we plot only the isochrones with ages between 
log(f/yr) = 7.0 and log(i/yr) = 10.5, in steps of Alogt = 
0.5. As already noted by Salaris & Weiss (1998), for solar 
and super-solar total mctallicity, the clement abundance 



ratios among metals affect the evolution. From the HR 
diagram of Fig. ||, two main features can be singled out. 

(i) cv-Enhanced isochrones have fainter and hotter 
turn-offs (TO). This is shown by the entries of Table |[ 
which lists the turn-off luminosity and effective tempera- 
ture of isochrones with Z = 0.019 but different mixtures of 
a-elements. The reason is that at given central hydrogen 
content the stellar model with solar-scaled composition 
has a mean opacity higher than the a-enhanced one. The 
point is illustrated in Fig. |^ where we compare the opac- 
ity profile across the 1 M Q models with [Z = 0.040, Y = 
0.320] and both a-enhanced (dashed line) and solar-scaled 
(solid line) mixtures. Higher opacities induce steeper ra- 
diative temperature gradients V r and lower surface tem- 
peratures T e g for the same central conditions. Further- 
more, steeper V r s imply smaller burning regions and lower 
luminosities in turn. This effect overwhelms the one in- 
duced by the slightly larger convective core due to the 
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Table 4. Surface chemical composition (by mass fraction) of the [Z — 0.008, Y — 0.250] solar-scaled models. 
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higher V r (when a convective core is present). Combining 
all those effects together, at the same evolutionary stage, 
the solar-scaled model is fainter, cooler and older than 
the a-enhanced one. The same trend is recovered in the 
isochrones as well. 

The implications of the different turn-off temperatures 
and luminosities on the ages of globular clusters have 
already been extensively investigated (see e.g. Salaris & 
Weiss 1998). For isochrones younger than 10 9 yr, the dif- 
ferences in the turn-off properties are much less marked. 
This is due both to (1) the weak metallicity dependence 
of the opacities in massive stars (where electron-scattering 
becomes the main source of opacity) and (2) the larger ex- 
tension of central convection caused by the flatter profile 
of V r in stellar interiors. 

(ii) a-Enhanced isochrones are hotter than the solar- 
scaled ones throughout all evolutionary phases. The ef- 
fect increases with the metallicity. It is mainly due to 
the higher opacities of the solar-scaled mixtures. Look- 
ing at the 10 Gyr isochrone with chemical composition 
[Z = 0.040, Y = 0.320] as an example, the temperature 



difference logT off solar - logT off Q ~ cnh is about -0.015 at 
the turn-off and increases to —0.020 at the bottom of the 
RGB. These differences in effective temperatures result 
into colour differences (B — V) so \ al — (B — V) a - cn h of 
0.05 mag and 0.08 mag, at the turn-off and bottom of 
the RGB, respectively. 

It is worth remarking that in Weiss et al. (1995), the a- 
enhanced isochrones are found to have hotter and brighter 
turn-offs, and cooler RGBs (at least in large portions of 
it) than in the solar-scaled case (see their Fig. 4; notice 
that the line-types there are wrong: solid lines refer to the 
a-enhanced case). The cooler RGBs are no longer found 
here or in VandenBerg et al. (2000) and have been a conse- 
quence of the lack of low-temperature opacity tables for a- 
enhanced compositions in opacities in Weiss et al. (1995). 
This confirms how important is to use self-consistent opac- 
ities. 

Finally, one can ask whether an a-enhanced isochrone 
may be approximated by a solar-scaled one with a slightly 
lower metallicity. In order to answer this question let us 
examine a set of 10 Gyr old isochrones. In Fig. |l^ we plot 
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log T e„ 

Fig. 8. Comparison of isochrones with chemical com- 
position [Z = 0.019, Y = 0.273] and different en- 
hancement of a-elements. Dotted lines refer to the a- 
enhanced isochrones, continuous lines to the solar-scaled 
ones. Only the isochrones with ages between log(i/yr) = 
7.0 and log(£/yr) = 10.5, in steps of Alogi = 0.5 for the 
sake of clarity, are plotted. 



all the isochrones with this age, and with metallicity going 
from Z = 0.008 to Z — 0.07 for both the a-enhanced (con- 
tinuous lines) and solar-scaled (short dashed lines) mix- 
tures. One can see that both sets of isochrones have dif- 
ferent shapes, the differences becoming larger at higher 
metallicities. More specifically, one can notice that the 
difference in effective temperatures at the turn-off point 
between an a-enhanced isochrone and a solar-scaled one, 
ATeff(TO), increases with the metallicity, passing from 
ATcfF(TO) = 0.008 to 0.015 as Z increases from 0.008 to 
0.07. On the contrary, the temperature difference at the 
base of the RGB, AT c g (RGB), is almost constant and of 
about AT eff (RGB) - 0.019. 

For the same age we calculate, by interpolating among 
the grids of different metallicity, several solar-scaled 
isochrones with Z spanning from 0.008 to 0.019, and 
selected among them the one best fitting the main se- 
quence of the Z = 0.019 a-enhanced isochrone. This 
turns out to be the interpolated isochrone with Z = 
0.011, which is shown by the long-dashed line in Fig. |l^. 
Clearly, this solar-scaled Z = 0.011 isochrone does not 
reproduce well the RGB sequence of the Z = 0.019 a- 
enhanced one. More specifically, although AT c g(TO) is 
very small for this isochrone pair, AT g (RGB) is not ne- 
glegible (AT eff (RGB) - 0.008). 




m/M 

Fig. 9. Comparison of the internal structure of two 1 M Q 
models with [Z — 0.040, Y = 0.320] at the same evolution- 
ary stage, i.e. when the central hydrogen mass fraction is 
0.36. The left scale refers both to the hydrogen abundance 
by mass and the opacity. The right scale refers to the log- 
arithm of the temperature. Continuous lines are for the 
solar-scaled model, whereas the dashed lines are for the 
a-enhanced one. 

From this simple test we then conclude that for 
a fixed age and relatively high metallicities, an a- 
enhanced isochrone cannot be reproduced by tuning the 
metallicity of a solar-scaled one. 

The above temperature differences for RGB models 
(isochrones) bear very much on integrated spectral line 
indices of old metal-rich populations, which are commonly 
used to infer ages and metallicities of elliptical galaxies. 
This will be the subject of a forthcoming paper. 

3.2. Theoretical Johnson and HST photometry 

Theoretical luminosities and effective temperatures are 
translated to magnitudes and colours by means of ex- 
tensive tabulations of bolometric corrections (BC) and 
colours obtained from properly convolving spectral energy 
distributions (SEDs) as a function of T e g , log g and [Fe/H] . 
The procedure is amply described in Bertelli et al. (1994), 
Bressan et al. (1994), and Tantalo et al. (1996) to whom 
the reader should refer for more details. Suffice it to recall 
here the following basic steps: 

(i) The main body of the spectral library is from Ku- 
rucz (1992), however extended to the high and low tem- 
perature ranges. For stars with T c g > 50, 000 K pure 
black-body spectra are assigned, whereas for stars with 
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Table 6. Position of the turn-off in the HR diagram of 
isochrones with metallicity Z = 0.019. 



solar-scaled a-enhanced 
age I log(L/L ) I logT cff I log(L/L ) I logT eff 



7.00 


4.491 


4.424 


4.473 


4.428 


7.20 


4.109 


4.367 


4.084 


4.370 


7.40 


3.749 


4.309 


3.722 


4.312 


7.60 


3.418 


4.253 


3.380 


4.255 


7.80 


3.095 


4.197 


3.046 


4.198 


8.00 


2.762 


4.143 


2.728 


4.143 


8.20 


2.390 


4.087 


2.367 


4.088 


8.40 


2.090 


4.030 


2.028 


4.031 


8.60 


1.791 


3.976 


1.733 


3.978 


8.80 


1.492 


3.924 


1.436 


3.926 


9.00 


1.207 


3.874 


1.154 


3.877 


9.20 


0.853 


3.834 


0.839 


3.837 


9.40 


0.587 


3.808 


0.528 


3.812 


9.60 


0.383 


3.787 


0.377 


3.795 


9.80 


0.379 


3.776 


0.406 


3.785 


10.00 


0.240 


3.759 


0.232 


3.770 


10.20 


0.083 


3.740 


0.082 


3.753 



T e ff < 3500 K the catalogue of stellar fluxes by Fluks 
et al. (1994) is adopted. This catalogue includes 97 ob- 
served spectra for all M-spectral subtypes in the wave- 
length range 3800 < A < 9000, and synthetic photo- 
spheric spectra in the range 9900 < A < 12500. 

The scale of T c g in Fluks et al. (1994) is similar to that 
of Ridgway et al. (1980) for spectral types earlier than M4 
but deviates from it for later spectral types. Since Ridg- 
way's et al. (1980) scale does not go beyond the spectral 
type M6, no comparison for more advanced spectral types 
is possible. 

The problem is further complicated by possible effects 
of metallicity. The Ridgway scale of T c g is based on stars 
with solar metallicity (Z ~ 0.02) and empirical calibra- 
tions of the T c ff-scale for Z ^ 0.02 are not available. 

To cope with this difficulty, we have introduced the 
metallicity-Teff relation of Bessell et al. (1989, 1991) using 
the (V — K) colour as a temperature indicator. An inter- 
polation is made between the T e g of Bessell et al. (1989) 
and the (V — K) colours given by Fluks et al. (1994) for 
the spectral types from M0 to M10. 

(ii) At a given metallicity [Fe/H], surface gravity log g, 
and T e g, BCs are determined by convolving the corre- 
sponding stellar SED with the response functions of the 
various pass-bands. 

(iii) We recall that Kurucz does not provide atmo- 
spheric models for both solar-scaled and a-enhanced mix- 
tures. It is not therefore possible to estimate the impact (if 
any) of a-enhanced atmospheres on the isochrone colours 
and magnitudes. 




Fig. 10. Comparison between a-enhanced (continuous 
lines) and solar-scaled (short-dashed lines) isochrones of 
10 Gyr. In both cases, from left to right the metallicity val- 
ues are Z = 0.008, 0.019, 0.04, and 0.07. We also present a 
solar-scaled isochrone with Z = 0.011 (long-dashed line), 
that has been obtained by interpolation among the solar- 
scaled tracks. This latter should be compared to the a- 
cnhanced isochrone with Z = 0.019 (see text). 

Johnson-Cousins photometry. The response functions for 
the various pass-bands in which magnitudes and colours 
are generated are from the following sources: Buser & Ku- 
rucz (1978) for the UBV pass-bands, Bessell (1990) for the 
R and / Cousins pass-bands, and finally Bessell & Brett 
(1988) for the JHK pass-bands. 

Re-normalization of the colours obtained by convolv- 
ing the SEDs with pass-band, has been made by convolv- 
ing the SED of Vega and imposing that the computed 
colours strictly match the observed ones (Kurucz 1992). 

Finally, the zero point of the BCs is fixed by imposing 
that the BC for the Kurucz model of the Sun is —0.08. 

HST-WFPC2 photometry. The transformation to the 
WFPC2 photometric system (filters F170W, F218W, 
F255W, F300W, F336W, F439W, F450W, F555W, 
F606W, F702W, F814W, and F850LP) requires some ad- 
ditional explanations. 

The main parameters defining this photometric system 
are given in Table [?]. Columns (1) to (3) report the filter 
name, the mean wavelength A and the r.m.s. band width 
(BW), respectively. 

With the WFPC2 detector at least three kinds of mag- 
nitudes are commonly in use: ttist, tivega and ttiab (see 
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Filter 


A (A) 


BW (A) 


Z P 


msT — JTlVEGA 


Note 


(1) 


(2) 


(3) ' 


(4) 


(5) 


(6) 



F170W 


1747 


290 


-20.6926 


-0.4070 




F218W 


2189 


171 


-20.8737 


-0.2263 




F255W 


2587 


170 


-21.0940 


-0.0060 




F300W 


2942 


325 


-21.1336 


+0.0336 


Wide U 


F336W 


3341 


204 


-21.1876 


+0.0876 


U-John. 


F439W 


4300 


202 


-20.4325 


-0.6675 


B-John. 


F450W 


4519 


404 


-20.6339 


-0.4661 


Wide B 


F555W 


5397 


522 


-21.0798 


-0.0202 


V-John. 


F606W 


5934 


637 


-21.4250 


+0.3250 


Wide V 


F702W 


6862 


587 


-21.8594 


+0.7594 


Wide R 


F814W 


7924 


647 


-22.3309 


+1.2309 


I- John. 


F850LP 


9070 


434 


-22.7578 


+1.6578 





the Synphot User's Guide by White et al. 1998, distributed 
by the STSDAS Group). In the following we limit our- 
selves to consider only the theoretical counterparts of tost 
and tovega- 

For radiation with flat spectrum and specific flux Fp 
(in ergs -1 cm -2 A ) impinging on the telescope, the 
tost magnitude is defined by 



TO3T,p = -21.1 - 2.5 log F P 



(1) 



which means that tost,p = for Fp = 3.63 10 9 
erg s" 1 cm -2 A . 

Let us then define the pass-band F(A) as the product 
of the filter transmission by the response function of the 
telescope assembly and detector in use. In this case, Fp is 
related to the counts rate Cp - i.e. the number of photons 
per second registered by the detector - through 



F F 



he 



A J P(X)XdX 



(2) 



where A is the effective collecting surface of the telescope, 
and h and c are the Planck constant and speed of light, 
respectively. 

Cp, the quantity actually measured at the telescope, 
can be easily converted to a specific flux by means of 
Eq. (||). Therefore, specific fluxes Fp and tost,p magni- 
tudes can be defined for incoming spectra of any shape. 
In the case of a stellar source of radius R, located at a 
distance d, and emitting a specific flux Ik, Eq. (|l|) can be 
replaced by 



tost,p = —21.1 — 2.5 log 



i? 2 fP(A)/ K (A)AdA 



d 2 J P(A) AdA 



(3) 



which defines the apparent tost.p magnitude of a star. 

In our case, we derive the absolute Mst.p magnitudes 
by simply locating our synthetic stars at a distance of 
d = 10 pc in Eq. The specific fluxes /k(A) are taken 



from the library of synthetic spectra in use as a function 
of T e ff, surface gravity g, and chemical composition. 
Similarly, we can define the tovega.p magnitudes 



tovega,p = Zp- 2.5 log 



R 2 J P(X)I K (X)\d\ 



d 2 J P(X) AdA 



(4) 



by imposing that the zero points Zp are such that the 
synthetic magnitudes of Vega match the ground-based 
UBVRI Johnson apparent magnitudes for this star, for 
the filters which are closest in wavelength to each other 
(see column 6 in Table [?]) . 

We adopt the values 0.02, 0.02, 0.03, 0.039, 0.035 mag 
for the apparent magnitudes of Vega in the U BVRI John- 
son system (see Holtzman et al. 1995). For the UV filters, 
Vega is assumed to have an apparent magnitude equal to 
zero. 

We remind the reader that if we assume the apparent 
Vega magnitude to be equal to zero in all the pass-bands, 
we would obtain Zp = 2.51ogFpv e ga and Eq. (Q) would 
become 



to V ega,p = -2.5 log 



Fp 



Fp,Veg 



(5) 



This latter is the definition of vegamag in the SYNPHOT 
package distributed with the STSDAS software. In this 
paper we do not make use of this calibration, but prefer 
to follow the slightly different one defined by Holtzman et 
al. (1995) and described above. 

To get the final calibration, we insert in Eq. ^ the spe- 
cific flux of Vega calculated by Castelli & Kurucz (1994) 
with atmosphere models, and converted it to an absolute 
flux on the Earth assuming for Vega the angular diam- 
eter of 3.24 mas (Code et al. 1976). The absolute flux 
we have adopted differs (a few percent) from the one by 
Hayes (1985) used by Holtzman et al. (1995). Once more, 
the absolute magnitude Mvega,p is obtained locating the 
synthetic star at a 10 pc distance. Column (4) of Table [?] 
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lists the zero points we have obtained, whereas column 
(5) gives the difference —21.1 — Zp between the zero point 
of the msT and ibvega magnitudes. The conversion from 
one magnitude scale to another is thus possible. 

Finally, we calculate the bolometric corrections BC 
which allow one to convert theoretical bolometric mag- 
nitudes both to Mst and Mvega: 

^(lOpc) 2 ^- 



M 



2.5 log 



21.1 + 2.51og 



BC(P) 



M ST 



L 

/P(A)J K (A)AdA 

/P(A)AdA 
21.1 - Z P 



(6) 
(7) 



BC(P)^f™ GA 

where a is the Stefan-Boltzmann constant. For the solar 
luminosity and absolute bolometric magnitude we adopt 



. - 3.844 10 33 erg s" 1 (Bahcall et al. 1995) and M° ol - 
4.77. 

Owing to the presence of contaminants inside the 
WFPC2 (see Baggett & Gonzaga 1998 and Holtzman et 
al. 1995), -P(A) changes slowly with time; in particular the 
UV throughput degrades, changing the photometric per- 
formances. To cope with this drawback, small corrections 
are usually added to the definition of the instrumental 
magnitudes in order to bring the magnitudes back to the 
optimal conditions (see Holtzman et al. 1995). Because of 
this, in the calculation of the bolometric corrections we 
do not use the present-day pass-bands P(A) provided by 
SYNPHOT but insert the pre-launch pass-bands, as in 
Holtzman et al. (1995). 

3.3. Integrated magnitudes and mass-to-light ratios of 
single stellar populations 

In this section we present the integrated magnitudes, 
colours and mass-to-light ratios for single stellar popu- 
lations (SSP) of the same age and chemical compositions 
of the isochrones above, both in the Johnson-Cousins and 
WFPC2 photometric systems. 

The integrated magnitudes and colours are computed 
assuming an initial mass function (f>(M) and adopting the 
normalization 

<-M„ 

M 4>{M) dM = 1 M (8) 

I Mi 

where M\ and M u are the lower and upper limits of zero 
age main sequence stars. By doing so, the integrated mag- 
nitudes of these ideal SSPs can be easily scaled to produce 
integrated magnitudes for stellar populations of any arbi- 
trary total initial mass Mx,i, by simply adding the factor 
-2.51ogM T ,i. 

The theoretical mass-to-light ratios M/L are calcu- 
lated following the same procedure as in Alongi & Chiosi 
(1990) and Chandar et al. (1999). 

The fate of a single star depends on its initial mass 
(binary stars are neglected here). In a brief and over- 
simplified picture of stellar evolution the following mass 



limits and groupings can be identified: (i) Stars more mas- 
sive than Mh ~ 0.7 — 0.8 M have a lifetime shorter than 
the current estimate of the Hubble age of the Universe, 
say 13 — 15 Gyr. Stars lighter than this limit are not of 
interest here because once formed they live for ever, (ii) 
Stars more massive than M up explode as supernovae leav- 
ing a neutron star remnant of 1.4 M®. The possibility that 
massive stars (M > 25 M ) may end up as Black Holes is 
not considered here, (iii) Stars less massive than Mw and 
more massive than Mh terminate their evolution as White 
Dwarfs of suitable masses M^d that depend on the initial 
mass and efficiency of mass loss during the RGB and AGB 
phases. With the current estimates of the mass loss effi- 
ciency Mw ~ 5 — 6M Q , i.e. Mw ~ M up . The possibility, 
however, exists that Mw < M up . In such a case stars in 
the mass interval Mw to M up reach the C-ignition stage 
and deflagrate as supernovae leaving no remnant. For the 
purposes of the present study, this latter case is neglected 
and Mw = M up is always assumed. Finally, let M a be the 
initial mass of the most massive star still alive in a SSP of 
a certain age. It is worth recalling that all the above mass 
limits depend on the initial chemical composition and that 
this dependence can be properly taken into account. 

Therefore, at any given age t the total mass Mr of a 
SSP is made of two contributions: 

pM a pM u 

M T (t)= / M (j>(M) dM + / M r 4>(M) dM (9) 

J Mi JM^ 

where the first right-hand side term is the total mass of 
the stars still alive at the age t (i.e. burning a nuclear fuel) 
and the second one is the total mass in stellar remnants, 
i.e. white dwarfs or neutron stars depending on the ini- 
tial mass of the progenitor (i.e. M r (M) = Mwd(-M) for 
M < M up , and M r (M) = 1.4 M for M > M up ). We 
adopt an upper mass limit M u = 120 M . The choice for 
the lower mass limit, Mi, is explained below. It is worth 
recalling here that, due to eq. ||, M-r(i) is the current mass 
of an ideal SSP with initial mass equal to 1 M . Mp(£) 
decreases as a function of time from its initial value, due 
to the mass lost by both stellar winds (for all masses) and 
supernova explosions (from massive stars). 

Finally, the total luminosity Lj of the SSP in any pass- 
band j is obtained by integrating, along the isochrone, the 
luminosity of individual stars whose number is given by 
the IMF 



M a 



Mi 



L j (M)<p{M) dM 



(10) 



To proceed further, one has to specify the initial mass 
function and M\. Here, we adopt the Salpeter (1955) law 
with slope x = 1.35 over the whole mass range M\ to M u . 
Concerning M\, the choice M\ = 0.068 M is made by 
imposing that the observational value MjLy = 0.19±0.04 
of the LMC cluster NGC 1866 is matched (see Girardi & 
Bica 1993). The age and metallicity of NGC 1866 are 10 8 
years and Z = 0.008, respectively, whereas its MjLy is 
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close to the mean value of 0.20 for LMC clusters of the 
same age (see Battinelli & Capuzzo-Dolcetta 1989, and 
references therein). 

4. Concluding remarks 

We have computed extended sets of a-enhanccd evolution- 
ary tracks and isochrones at relatively high mctallicities. A 
major improvement, with respect to previous calculations 
of a-enhanced tracks, is that we made use of self-consistent 
opacities, i.e. computed with the same a-enhanced chemi- 
cal compositions over the complete range of temperatures. 
The main result is that in general all evolutionary phases 
(isochrones) have higher effective temperature with re- 
spect to solar-scaled models (isochrones) of same metal 
content Z. The temperature shift is caused by the lower 
opacities of the a-enhanced mixtures. Moreover, we find 
that at relatively high mctallicities and old ages, an a- 
cnhanced isochrone cannot be mimicked by simply using 
a solar-scaled isochrone of lower metallicity. 

Obvious limitation of the present results is that they 
refer to a particular choice for the a-enhanced chemi- 
cal compositions derived from observations of a particu- 
lar sample of metal-poor field stars. Future observational 
data may suggest a different partition of metals in a- 
enhanced stars, and hence different results for the stellar 
models and corresponding isochrones. 

Testing a large number of possible combinations of a- 
enhanced ratios, though feasible, would not be of practical 
use owing to the present uncertainty on the abundance ra- 
tios. Therefore, the results of this study ought to be taken 
as indicative of the potential effect of a-enhanced element 
ratios. 

Moreover, it is likely that variations of single a- 
elements are less important than the whole problem of 
including or not the enhancement of a-elements in the 
initial chemical composition. Finally, it is still not clear 
whether differences in enhancements are significant com- 
pared to the errors in the determinations. 

Retrieval of the data sets. Complete tabulations of the 
relevant data for all evolutionary sequences, isochrones, 
integrated magnitudes and colours and mass-to-light ra- 
tions together with useful summary tables can be obtained 
either upon request to the authors or downloaded from the 
address http : //pie iadi .pd. astro . it. 

The layout of the tables of stellar models and 
isochrones is the same as in Girardi et al. (2000) as far 
as the Johnson-Cousins photometric system is concerned. 
For the WFPC2 system the layout is similar, with tabic 
headers allowing to easily identify the pass-band. 
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